#Setting a reasonable p-value threshold to be used throughout
p_cutoff <- 0.05
p_cutoff_new <- 0.05
knitr::opts_chunk$set(echo = FALSE, message = FALSE, warning = FALSE)
FC_cutoff_original <- 0
FC_cutoff_new <- 1
Global parameters for modeling
Exploratory Plots
Heatmaps
Individuals clustered by overall expression

Individuals by Top 100 genes heatmap

Individuals clustered by top 100 gene expression

PCA for overall expression

LS0tCnRpdGxlOiAiREVhbmFseXNpc19rZmlzaF9vc21vdGljIgphdXRob3I6ICJMaXNhIEpvaG5zb24iCmRhdGU6ICdgciBTeXMuRGF0ZSgpYCcKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICBjb2RlX2ZvbGRpbmc6IGhpZGUKICAgIGNvbGxhcHNlZDogbm8KICAgIGRmX3ByaW50OiBwYWdlZAogICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMKICAgIHRoZW1lOiBjZXJ1bGVhbgogICAgdG9jOiB5ZXMKICAgIHRvY19kZXB0aDogNQogICAgdG9jX2Zsb2F0OiB5ZXMKICBodG1sX25vdGVib29rOgogICAgdG9jOiB5ZXMKICAgIHRvY19kZXB0aDogNQotLS0KCmBgYHtyIEdsb2JhbFZhcmlhYmxlc30KCiNTZXR0aW5nIGEgcmVhc29uYWJsZSBwLXZhbHVlIHRocmVzaG9sZCB0byBiZSB1c2VkIHRocm91Z2hvdXQKCnBfY3V0b2ZmIDwtIDAuMDUKCnBfY3V0b2ZmX25ldyA8LSAwLjA1Cgprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IEZBTFNFLCBtZXNzYWdlID0gRkFMU0UsIHdhcm5pbmcgPSBGQUxTRSkKCkZDX2N1dG9mZl9vcmlnaW5hbCA8LSAwCgpGQ19jdXRvZmZfbmV3IDwtIDEKCgpgYGAKCgojR2xvYmFsIHBhcmFtZXRlcnMgZm9yIG1vZGVsaW5nCgoKCmBgYHtyIExvYWRQYWNrYWdlcywgcmVzdWx0cz0naGlkZScsIGluY2x1ZGU9RkFMU0V9CgojIEluc3RhbGwgZnVuY3Rpb24gZm9yIHBhY2thZ2VzICAgIApwYWNrYWdlczwtZnVuY3Rpb24oeCl7CiAgeDwtYXMuY2hhcmFjdGVyKG1hdGNoLmNhbGwoKVtbMl1dKQogIGlmICghcmVxdWlyZSh4LGNoYXJhY3Rlci5vbmx5PVRSVUUpKXsKICAgIGluc3RhbGwucGFja2FnZXMocGtncz14LHJlcG9zPSJodHRwOi8vY3Jhbi5yLXByb2plY3Qub3JnIikKICAgIHJlcXVpcmUoeCxjaGFyYWN0ZXIub25seT1UUlVFKQogIH0KfQoKYmlvY29uZHVjdG9ycyA8LSBmdW5jdGlvbih4KXsKICAgIHg8LSBhcy5jaGFyYWN0ZXIobWF0Y2guY2FsbCgpW1syXV0pCiAgICBpZiAoIXJlcXVpcmUoeCwgY2hhcmFjdGVyLm9ubHkgPSBUUlVFKSl7CiAgICAgIHNvdXJjZSgiaHR0cHM6Ly9iaW9jb25kdWN0b3Iub3JnL2Jpb2NMaXRlLlIiKQogICAgICBiaW9jTGl0ZShwa2dzPXgpCiAgICAgIHJlcXVpcmUoeCwgY2hhcmFjdGVyLm9ubHkgPSBUUlVFKQogICAgfQp9CgpwYWNrYWdlcyhNQVNTKQpwYWNrYWdlcyhnZ3Bsb3QyKQpwYWNrYWdlcyhndG9vbHMpCnBhY2thZ2VzKHBoZWF0bWFwKQpwYWNrYWdlcyhjb3dwbG90KQpwYWNrYWdlcyhSQ29sb3JCcmV3ZXIpCnBhY2thZ2VzKGRwbHlyKQpwYWNrYWdlcyh0aWR5cikKcGFja2FnZXMoZ2dyZXBlbCkKYmlvY29uZHVjdG9ycyhERVNlcTIpCmJpb2NvbmR1Y3RvcnMobGltbWEpCmJpb2NvbmR1Y3RvcnMoJ2VkZ2VSJykKcGFja2FnZXMoZ3Bsb3RzKQpwYWNrYWdlcyhsYXR0aWNlKQoKc2Vzc2lvbkluZm8oKQoKYGBgCgoKYGBge3IgbG9hZGZpbGVzLCByZXN1bHRzPSdoaWRlJywgaW5jbHVkZT1GQUxTRX0KCiMgVGhpcyBpcyB0aGUgY291bnRzIHdpdGggRXhwZXJpbWVudGFsIERlc2lnbiBJbmZvIGluIHRoZSBsYXN0IDUgcm93cwoKaWYoIWZpbGUuZXhpc3RzKCcuLi8uLi8uLi9FbnNlbWJsX3NwZWNpZXNfY291bnRzX2Rlc2lnbmZhY3RvcnMuY3N2JykpewogIGRvd25sb2FkLmZpbGUoImh0dHBzOi8vb3NmLmlvLzd2cDM4L2Rvd25sb2FkIiwnRW5zZW1ibF9zcGVjaWVzX2NvdW50c19kZXNpZ25mYWN0b3JzLmNzdicpCn0KaWYoIWZpbGUuZXhpc3RzKCcuLi8uLi8uLi9ub256ZXJvX2NsYWRlX3BoeXNpb2xvZ3lfY291bnRzX2Rlc2lnbi5jc3YnKSl7CiAgZG93bmxvYWQuZmlsZSgiaHR0cHM6Ly9vc2YuaW8vYmU3bnkvZG93bmxvYWQiLCdub256ZXJvX2NsYWRlX3BoeXNpb2xvZ3lfY291bnRzX2Rlc2lnbi5jc3YnKQp9CmlmKCFmaWxlLmV4aXN0cygnLi4vLi4vLi4vZ3JlYXRlcjVjb3VudHNfY2xhZGVfcGh5c2lvbG9neV9jb3VudHNfZGVzaWduLmNzdicpKXsKICBkb3dubG9hZC5maWxlKCJodHRwczovL29zZi5pby9iZTdueS9kb3dubG9hZCIsJ2dyZWF0ZXI1Y291bnRzX2NsYWRlX3BoeXNpb2xvZ3lfY291bnRzX2Rlc2lnbi5jc3YnKQp9CgojY291bnRzX2Rlc2lnbiA8LSByZWFkLmNzdigiRW5zZW1ibF9zcGVjaWVzX2NvdW50c19kZXNpZ25mYWN0b3JzLmNzdiIsc3RyaW5nc0FzRmFjdG9ycyA9IFRSVUUpCmNvdW50c19kZXNpZ24gPC0gcmVhZC5jc3YoIi4uLy4uLy4uL25vbnplcm9fY2xhZGVfcGh5c2lvbG9neV9jb3VudHNfZGVzaWduLmNzdiIsc3RyaW5nc0FzRmFjdG9ycyA9IFRSVUUpCgoKYGBgCgoKYGBge3IgZGVzaWduaW5mbywgaW5jbHVkZT1GQUxTRX0KIyAtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojIEZvcm1hdCBkZXNpZ24gYW5kIGNvdW50cyBtYXRyaXgKIyAtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKZGVzaWduIDwtIGNvdW50c19kZXNpZ25bY291bnRzX2Rlc2lnbiRFbnNlbWJsID09ICdFbXB0eScsXQojZGVzaWduJHR5cGUgPC0gYygic3BlY2llcyIsIm5hdGl2ZV9zYWxpbml0eSIsImNsYWRlIiwiZ3JvdXAiLCJjb25kaXRpb24iKQpkcm9wcyA8LSBjKCJYIiwiRW5zZW1ibCIsCiAgICAgICAgICAgIkZfemVicmludXNfQldfMS5xdWFudCIsIkZfemVicmludXNfQldfMi5xdWFudCIsCiAgICAgICAgICAgIkZfemVicmludXNfRldfMS5xdWFudCIsIkZfemVicmludXNfRldfMi5xdWFudCIsCiAgICAgICAgICAgIkZfbm90dGlfRldfMS5xdWFudCIsIkZfbm90dGlfRldfMi5xdWFudCIsCiAgICAgICAgICAgIkZfc2NpYWRpY3VzX0JXXzEucXVhbnQiLCJGX3NjaWFkaWN1c19GV18xLnF1YW50IiwiRl9zY2lhZGljdXNfRldfMi5xdWFudCIpCnRyYW5zZmVyX2Ryb3BzIDwtIGMoIkZfc2NpYWRpY3VzX3RyYW5zZmVyXzEucXVhbnQiLCJGX3JhdGhidW5pX3RyYW5zZmVyXzEucXVhbnQiLCJGX3JhdGhidW5pX3RyYW5zZmVyXzIucXVhbnQiLCJGX3JhdGhidW5pX3RyYW5zZmVyXzMucXVhbnQiLAogICAgICAgICAgICAgICAgICAgICJGX2dyYW5kaXNfdHJhbnNmZXJfMS5xdWFudCIsIkZfZ3JhbmRpc190cmFuc2Zlcl8yLnF1YW50IiwiRl9ncmFuZGlzX3RyYW5zZmVyXzMucXVhbnQiLAogICAgICAgICAgICAgICAgICAgICJGX25vdGF0dXNfdHJhbnNmZXJfMS5xdWFudCIsIkZfbm90YXR1c190cmFuc2Zlcl8yLnF1YW50IiwiRl9ub3RhdHVzX3RyYW5zZmVyXzMucXVhbnQiLAogICAgICAgICAgICAgICAgICAgICJGX3BhcnZhcGluaXNfdHJhbnNmZXJfMS5xdWFudCIsIkZfcGFydmFwaW5pc190cmFuc2Zlcl8yLnF1YW50IiwKICAgICAgICAgICAgICAgICAgICAiTF9nb29kZWlfdHJhbnNmZXJfMS5xdWFudCIsIkxfZ29vZGVpX3RyYW5zZmVyXzIucXVhbnQiLCJMX2dvb2RlaV90cmFuc2Zlcl8zLnF1YW50IiwKICAgICAgICAgICAgICAgICAgICAiRl9vbGl2YWNlb3VzX3RyYW5zZmVyXzEucXVhbnQiLCJGX29saXZhY2VvdXNfdHJhbnNmZXJfMi5xdWFudCIsCiAgICAgICAgICAgICAgICAgICAgIkxfcGFydmFfdHJhbnNmZXJfMS5xdWFudCIsIkxfcGFydmFfdHJhbnNmZXJfMi5xdWFudCIsIkxfcGFydmFfdHJhbnNmZXJfMy5xdWFudCIsCiAgICAgICAgICAgICAgICAgICAgIkZfaGV0ZXJvY2xpdHVzTURQUF90cmFuc2Zlcl8xLnF1YW50IiwiRl9oZXRlcm9jbGl0dXNNRFBQX3RyYW5zZmVyXzIucXVhbnQiLCJGX2hldGVyb2NsaXR1c01EUFBfdHJhbnNmZXJfMy5xdWFudCIsCiAgICAgICAgICAgICAgICAgICAgIkZfc2ltaWxpc190cmFuc2Zlcl8xLnF1YW50IiwiRl9zaW1pbGlzX3RyYW5zZmVyXzIucXVhbnQiLCJGX3NpbWlsaXNfdHJhbnNmZXJfMy5xdWFudCIsCiAgICAgICAgICAgICAgICAgICAgIkZfZGlhcGhhbnVzX3RyYW5zZmVyXzEucXVhbnQiLCJGX2RpYXBoYW51c190cmFuc2Zlcl8yLnF1YW50IiwKICAgICAgICAgICAgICAgICAgICAiRl9jaHJ5c290dXNfdHJhbnNmZXJfMS5xdWFudCIsIkZfY2hyeXNvdHVzX3RyYW5zZmVyXzIucXVhbnQiLAogICAgICAgICAgICAgICAgICAgICJBX3hlbmljYV90cmFuc2Zlcl8xLnF1YW50IiwiQV94ZW5pY2FfdHJhbnNmZXJfMi5xdWFudCIsIkFfeGVuaWNhX3RyYW5zZmVyXzMucXVhbnQiICwKICAgICAgICAgICAgICAgICAgICAiRl9jYXRhbmF0dXNfdHJhbnNmZXJfMS5xdWFudCIsIkZfY2F0YW5hdHVzX3RyYW5zZmVyXzIucXVhbnQiLAogICAgICAgICAgICAgICAgICAgICJGX2hldGVyb2NsaXR1c01EUExfdHJhbnNmZXJfMS5xdWFudCIsIkZfaGV0ZXJvY2xpdHVzTURQTF90cmFuc2Zlcl8yLnF1YW50IiwiRl9oZXRlcm9jbGl0dXNNRFBMX3RyYW5zZmVyXzMucXVhbnQiKQpjb3VudHM8LWNvdW50c19kZXNpZ25bIWNvdW50c19kZXNpZ24kRW5zZW1ibCA9PSAnRW1wdHknLF0Kcm93bmFtZXMoY291bnRzKTwtY291bnRzJEVuc2VtYmwKZGVzaWduIDwtIGRlc2lnblsgLCAhKG5hbWVzKGRlc2lnbikgJWluJSBkcm9wcyldCmNvdW50cyA8LSBjb3VudHNbICwgIShuYW1lcyhjb3VudHMpICVpbiUgZHJvcHMpXQpkZXNpZ24gPC0gZGVzaWduWyAsICEobmFtZXMoZGVzaWduKSAlaW4lIHRyYW5zZmVyX2Ryb3BzKV0KY291bnRzIDwtIGNvdW50c1sgLCAhKG5hbWVzKGNvdW50cykgJWluJSB0cmFuc2Zlcl9kcm9wcyldCmRpbShkZXNpZ24pCiNbMV0gIDUgODEKZGltKGNvdW50cykKI1sxXSAzMDQ2NiAgICA4MQojIC0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgZGVzaWduIGNhdGVvZ29yaWVzCiMgLS0tLS0tLS0tLS0tLS0tLS0tLS0KCnNwZWNpZXM8LWFzLmNoYXJhY3Rlcih1bmxpc3QoZGVzaWduWzEsXSkpCnBoeXNpb2xvZ3k8LWFzLmNoYXJhY3Rlcih1bmxpc3QoZGVzaWduWzIsXSkpCmNsYWRlPC1hcy5jaGFyYWN0ZXIodW5saXN0KGRlc2lnblszLF0pKQpjb25kaXRpb248LWFzLmNoYXJhY3Rlcih1bmxpc3QoZGVzaWduWzUsXSkpCmNvbmRpdGlvbl9waHlzaW9sb2d5PC1hcy52ZWN0b3IocGFzdGUoY29uZGl0aW9uLHBoeXNpb2xvZ3ksc2VwPSIuIikpCmNvbHM8LWNvbG5hbWVzKGNvdW50cykKRXhwRGVzaWduIDwtIGRhdGEuZnJhbWUocm93Lm5hbWVzPWNvbHMsCiAgICAgICAgICAgICAgICAgICAgICAgIGNvbmRpdGlvbj1jb25kaXRpb24sCiAgICAgICAgICAgICAgICAgICAgICAgIHBoeXNpb2xvZ3kgPSBwaHlzaW9sb2d5LAogICAgICAgICAgICAgICAgICAgICAgICBjbGFkZSA9IGNsYWRlLAogICAgICAgICAgICAgICAgICAgICAgICBzcGVjaWVzID0gc3BlY2llcywKICAgICAgICAgICAgICAgICAgICAgICAgc2FtcGxlPWNvbHMpCkV4cERlc2lnbgpkZXNpZ24gPSBtb2RlbC5tYXRyaXgoIH4gcGh5c2lvbG9neSArIGNvbmRpdGlvbiArIHBoeXNpb2xvZ3k6Y29uZGl0aW9uLCBFeHBEZXNpZ24pCgpjb2xuYW1lcyhkZXNpZ24pCiMgY2hlY2sgcmFuayBvZiBtYXRyaXgKTWF0cml4OjpyYW5rTWF0cml4KCBkZXNpZ24gKQpkaW0oZGVzaWduKQpgYGAKCgpgYGB7ciBub3JtLCByZXN1bHRzPSdoaWRlJywgaW5jbHVkZT1GQUxTRX0KY291bnRzX3JvdW5kIDwtIHJvdW5kKGRhdGEubWF0cml4KGNvdW50cyksIGRpZ2l0cz0wKQojY291bnRzX3JvdW5kIDwtIGhlYWQoY291bnRzX3JvdW5kLCBuID0gMTkwMDApCmRpbShjb3VudHNfcm91bmQpCiNwbG90KGNvbFN1bXModChjb3VudHNfcm91bmQpKSApCmxjcG0yIDwtIGNwbShjb3VudHNfcm91bmQsIGxvZyA9IFRSVUUpCiNib3hwbG90KGxjcG0yLCBsYXMgPSAyLCBtYWluID0gIkJlZm9yZSBERVNlcTIgTm9ybWFsaXphdGlvbiIpCgojIHRyYW5zZm9ybSB3aXRoIERFU2VxIHRvIHN0YWJpbGl6ZSB2YXJpYW5jZQoKZGRzIDwtIERFU2VxRGF0YVNldEZyb21NYXRyaXgoY291bnREYXRhID0gY291bnRzX3JvdW5kLGNvbERhdGEgPSBFeHBEZXNpZ24sZGVzaWduID0gZGVzaWduKQpkZHMgPC0gZXN0aW1hdGVTaXplRmFjdG9ycyhkZHMpCmRkcyA8LSBlc3RpbWF0ZURpc3BlcnNpb25zKGRkcykKZGRzJHNpemVGYWN0b3IKbmMgPC0gY291bnRzKGRkcywgbm9ybWFsaXplZD1UUlVFKQojd3JpdGUuY3N2KG5jLCAiLi4vLi4vLi4vbm9ybWNvdW50cy5jc3YiKQpjb3VudHM8LW5jCiNwbG90KGNvbFN1bXModChjb3VudHMpKSApCgpjcG0gPC0gY3BtKGNvdW50cykKbGNwbSA8LSBjcG0oY291bnRzLCBsb2cgPSBUUlVFKQojcGxvdChjb2xTdW1zKHQobGNwbSkpICkKI2JveHBsb3QobGNwbSwgbGFzID0gMiwgbWFpbiA9ICJBZnRlciBERVNlcTIgTm9ybWFsaXphdGlvbiIpCmBgYAoKCiMgRXhwbG9yYXRvcnkgUGxvdHMKCgojIyBIZWF0bWFwcwoKCiMjIyBJbmRpdmlkdWFscyBjbHVzdGVyZWQgYnkgb3ZlcmFsbCBleHByZXNzaW9uCgoKYGBge3IgUGxhaW5IZWF0bWFwLCBmaWcua2VlcD0ibGFzdCIsIGZpZy53aWR0aD0xMSwgZmlnLnBhdGg9J2ZpZ3VyZXMvJywgZGV2PWMoJ3BuZycsICdwZGYnKX0KCgoKcmxkIDwtIHZzdChkZHMsIGJsaW5kID0gRkFMU0UsZml0VHlwZT0nbG9jYWwnKQpzYW1wbGVEaXN0cyA8LSBkaXN0KHQoYXNzYXkocmxkKSkpCmRmIDwtIGFzLmRhdGEuZnJhbWUoY29sRGF0YShkZHMpWyxjKCJwaHlzaW9sb2d5IiwiY29uZGl0aW9uIiwiY2xhZGUiKV0pCnNhbXBsZURpc3RNYXRyaXggPC0gYXMubWF0cml4KCBzYW1wbGVEaXN0cyApCmNvbG9ycyA8LSBjb2xvclJhbXBQYWxldHRlKCByZXYoYnJld2VyLnBhbCg5LCAiQmx1ZXMiKSkgKSgyNTUpCnBoZWF0bWFwKHNhbXBsZURpc3RNYXRyaXgsCiAgICAgICAgIGNsdXN0ZXJpbmdfZGlzdGFuY2Vfcm93cyA9IHNhbXBsZURpc3RzLAogICAgICAgICBjbHVzdGVyaW5nX2Rpc3RhbmNlX2NvbHMgPSBzYW1wbGVEaXN0cywKICAgICAgICAgY29sID0gY29sb3JzLCBhbm5vdGF0aW9uID0gZGYsIHNob3dfcm93bmFtZXM9RikKCmBgYAoKIyMjIEluZGl2aWR1YWxzIGJ5IFRvcCAxMDAgZ2VuZXMgaGVhdG1hcAoKCgpgYGB7ciBNaW5pUGxhaW5HZW5lSGVhdG1hcCwgZWNobz1GQUxTRSwgZmlnLmtlZXA9Imxhc3QiLCBmaWcud2lkdGg9MTEsIGZpZy5wYXRoPSdmaWd1cmVzLycsIGRldj1jKCdwbmcnLCAncGRmJyl9CgpzZWxlY3QxMDAgPC0gb3JkZXIocm93TWVhbnMoY291bnRzKGRkcyxub3JtYWxpemVkPVRSVUUpKSxkZWNyZWFzaW5nPVRSVUUpWzE6MTAwXQoKcGhlYXRtYXAoYXNzYXkocmxkKVtzZWxlY3QxMDAsXSwgc2hvd19yb3duYW1lcz1ULGNsdXN0ZXJpbmdfZGlzdGFuY2Vfcm93cyA9IHNhbXBsZURpc3RzLAogICAgICAgICBjbHVzdGVyaW5nX2Rpc3RhbmNlX2NvbHMgPSBzYW1wbGVEaXN0cywgYW5ub3RhdGlvbl9jb2w9ZGYpCgoKYGBgCgojIyMgSW5kaXZpZHVhbHMgY2x1c3RlcmVkIGJ5IHRvcCAxMDAgZ2VuZSBleHByZXNzaW9uCgoKYGBge3IgTWluaVBsYWluSGVhdG1hcCwgZWNobz1GQUxTRSwgZmlnLmtlZXA9Imxhc3QiLCBmaWcud2lkdGg9MTEsIGZpZy5wYXRoPSdmaWd1cmVzLycsIGRldj1jKCdwbmcnLCAncGRmJyl9CgpzZWxlY3QxMDAgPC0gb3JkZXIocm93TWVhbnMoY291bnRzKGRkcyxub3JtYWxpemVkPVRSVUUpKSxkZWNyZWFzaW5nPVRSVUUpWzE6MTAwXQoKc2FtcGxlRGlzdHMgPC0gZGlzdCh0KGFzc2F5KHJsZClbc2VsZWN0MTAwLF0pKQpzYW1wbGVEaXN0TWF0cml4IDwtIGFzLm1hdHJpeCggc2FtcGxlRGlzdHMgKQoKcGhlYXRtYXAoc2FtcGxlRGlzdE1hdHJpeCwgc2hvd19yb3duYW1lcz1ULGNsdXN0ZXJpbmdfZGlzdGFuY2Vfcm93cyA9IHNhbXBsZURpc3RzLAogICAgICAgICBjbHVzdGVyaW5nX2Rpc3RhbmNlX2NvbHMgPSBzYW1wbGVEaXN0cywgYW5ub3RhdGlvbl9jb2w9ZGYpCgoKYGBgCgojIyMgUENBIGZvciBvdmVyYWxsIGV4cHJlc3Npb24KCgpgYGB7ciBwbGFpblBDQSwgZmlnLmtlZXA9Imxhc3QiLCBmaWcud2lkdGg9MTEsIGZpZy5wYXRoPSdmaWd1cmVzLycsIGRldj1jKCdwbmcnLCAncGRmJyl9CgoKY293cGxvdDo6cGxvdF9ncmlkKCBwbG90UENBKHJsZCwgaW50Z3JvdXA9ImNvbmRpdGlvbiIpLAogICAgICAgICAgICAgICAgICAgIHBsb3RQQ0EocmxkLCBpbnRncm91cD0icGh5c2lvbG9neSIpLAogICAgICAgICAgICAgICAgICAgIHBsb3RQQ0EocmxkLCBpbnRncm91cD0iY2xhZGUiKSwKICAgICAgICAgICAgICAgICAgICBwbG90UENBKHJsZCwgaW50Z3JvdXA9YygiY2xhZGUiLCJwaHlzaW9sb2d5IikpLAogICAgICAgICAgICAgICAgICAgICAgICAgICBhbGlnbj0iYyIsIG5jb2w9MikKCmBgYAoK